Perturbation theory for self-gravitating gauge fields I: The odd-parity sector 
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A gauge- and coordinate-invariant perturbation theory for self-gravitating non-Abelian gauge fields 
is developed and used to analyze local uniqueness and linear stability properties of non-Abelian 
equilibrium configurations. It is shown that all admissible stationary odd-parity excitations of the 
static and spherically symmetric Einstein- Yang-Mills soliton and black hole solutions have total 
angular momentum number 1=1, and are characterized by non- vanishing asymptotic flux inte- 
grals. Local uniqueness results with respect to non-Abelian perturbations are also established for 
the Schwarzschild and the Reissner-Nordstrom solutions, which, in addition, are shown to be linearly 
stable under dynamical Einstein- Yang-Mills perturbations. Finally, unstable modes with 1=1 are 
also excluded for the static and spherically symmetric non-Abelian solitons and black holes. (PACS 
numbers: 04.25.Nx, 04.40.-b, 04.70.Bw) 



I. INTRODUCTION 

Self-gravitating non-Abelian gauge fields admit a rich 
spectrum of equilibrium configurations, which is a conse- 
quence of the balance between the gravitational attrac- 
tion and the repulsive nature of the Yang-Mills interac- 
tion. In particular, the static and spherically symmetric 
non-Abelian soliton [Q and black hole solutions owe 
their existence to the nonlinearities of both general rela- 
tivity and Yang-Mills theory. 

On the other hand, the key to the black hole unique- 
ness theorems Q lies in the cr-model structure of the 
Einstein(-Maxwell) equations in the presence of a Killing 
field [Q, 0. As this property ceases to exist for self- 
gravitating non-Abelian gauge fields [pj, the classification 
of all stationary Einstein- Yang-Mills (EYM) soliton and 
black hole solutions is necessarily a very difficult task. In 
particular, the set of global charges (asymptotic flux in- 
tegrals) does no longer uniquely characterize all possible 
non-Abelian equilibrium configurations. 

Induced by the work of Bartnik and McKinnon (BK) 
on non-Abelian solitons m, various new self-gravitating 
equilibrium configurations have been found during the 
last decade. Besides the abovementioned static and 
spherically symmetric black holes with Yang-Mills hair 
(i.e., with vanishing Yang-Mills charges but different met- 
ric structure than the Schwarzschild solution) |Q], these 
include soliton and black hole solutions in Skyrme, Higgs, 
dilaton and other non-linear field theories coupled to 
gravity (see M for a review and references). 

Moreover, numerical ||] and analytical Q studies have 
revealed that non-Abelian static black holes are not nec- 
essarily spherically symmetric - in fact, they need not 
even be axisymmetric [ |lO| . In addition, the non- linear 
nature of the Yang-Mills interaction enables the exis- 
tence of stationary, non-static black holes with vanishing 
Komar angular momentum |11| . Also, the usual Lewis- 
Papapetrou form of the metric does not necessarily de- 
scribe all stationary and axisymmetric EYM black holes, 
that is, the circularity theorem does not generalize to 



space-times containing non-Abelian gauge fields [Of. 

The above comments suggest that it is not (yet) fea- 
sible to completely classify the soliton and black hole 
solutions of the stationary EYM equations. In this ar- 
ticle we pursue, therefore, a more modest aim. That is, 
we compute the complete spectrum of stationary EYM 
perturbations of the BK solitons and the corresponding 
black holes with hair. We do so by systematically devel- 
oping the perturbation theory for self-gravitating non- 
Abelian gauge fields. Following the tradition, we start 
with the odd-parity sector, and defer the investigation of 
even-parity perturbations to a forthcoming publication 

The gauge- and coordinate-invariant equations derived 
in this paper describe perturbations of arbitrary spheri- 
cally symmetric EYM configurations, where the station- 
ary and the dynamical sector can be treated separately 
if the background is static. In order to classify the equi- 
librium solutions close to the BK solitons and the corre- 
sponding black holes, it is sufficient to consider stationary 
excitations. As we shall see, these are naturally analyzed 
in terms of invariant metric and Yang-Mills amplitudes. 

The main results of this paper concern two local 
uniqueness theorems, applying to the BK solitons and 
the corresponding black holes with hair, respectively: 
We prove that all stationary odd-parity excitations of 
these static and spherically symmetric configurations are 
parametrized in terms of infinitesimal asymptotic flux 
integrals. More precisely, we show that the soliton and 
black hole excitations found in |ll| are the only station- 
ary, asymptotically flat perturbations of the BK solitons 
and the corresponding black holes with hair. In partic- 
ular, there exist no admissible regular or black hole per- 
turbations with total angular momentum number £ > 1, 
while for i = 1, the unique soliton and black hole excita- 
tions are those with infinitesimal electric charge and/or 
infinitesimal Komar angular momentum 1 11 1. On the per- 
turbative level, the situation is, therefore, similar to the 
Abelian case, where the only admissible stationary exci- 
tations of the Schwarzschild metric are the Kerr-Newman 
modes. The above results also establish a local version 



of the circularity theorem in the odd-parity sector. 

In addition to the classification of neighboring equi- 
librium configurations, we also discuss some stability is- 
sues, which require the analysis of dynamical perturba- 
tions. Unfortunately, the gauge-invariant metric pertur- 
bations used in this paper are, in general, not suited 
to apply spectral analysis, since their evolution is not 
governed by a standard pulsation operator. In a recent 
work |14] we have demonstrated how to overcome this 
problem by using curvature-based quantities. A rigorous 
discussion of dynamical perturbations within the metric 
approach is nevertheless possible for some distinguished 
cases. These include i = 1 EYM perturbations of ar- 
bitrary background configuratiuons, and arbitrary EYM 
perturbations of embedded Abelian configurations. 

Hence, further results derived in this paper concern 
the non- Abelian stability (and local uniqueness) of the 
Schwarzschild and the Reissner-Nordstrom (RN) black 
holes, as well as the stability properties of non- Abelian 
configurations with respect to £ = 1 perturbations. In 
particular, we show that both the Schwarzschild and the 
RN metric are linearly stable with respect to dynamical 
non- Abelian perturbations and admit no stationary exci- 
tations other than the (embedded) Kerr-Newman modes. 
In addition, we establish the absence of unstable modes 
of the pulsation equations governing the £ = 1 pertur- 
bations of the BK solitons and the corresponding black 
holes with hair. In this context it is worthwhile recalling 
that unstable Yang-Mills modes with odd-parity do exist 
for^ = |l|. 

The paper is organized as follows: In Section f| 
we briefly review the gauge-invariant approach to odd- 
parity gravitational perturbations and give a coordinate- 
invariant derivation of the Regge- Wheeler (RW) equa- 
tion. In Section p| we present the harmonic decompo- 
sition of Yang-Mills fields, using a convenient method to 
parametrize su(2)-valued one-forms in terms of isospin 
harmonics. Taking advantage of some powerful tools de- 
veloped in Appendix y|, the linearized field equations 
governing arbitrary odd-parity perturbations of spheri- 
cally symmetric EYM configurations are derived in Sec- 
tion 



[V 



As first applications, we establish the linear stability 
and the local uniqueness properties of the Schwarzschild 
and the RN solutions with respect to non-Abelian per- 
turbations in Sections M and ^, respectively. The local 
uniqueness theorems for the BK solitons a nd t he corre- 
sponding black holes are proven in Section VII. Eventu- 



ally, in Section VIII , we establish the dynamical stability 
of these solutions with respect to non-spherical pertur- 
bations with £ — 1. 

A variety of technical issues, such as the expressions 
for the linearized Ricci tensor, the integral argument ex- 
cluding admissible solutions of certain RW type equa- 
tions, some asymptotic expansions, the introduction of 
isospin harmonics, and the construction of gauge- and 
coordinate-invariant Yang-Mills amplitudes are discussed 
in Appendixes \ 



II. GRAVITATIONAL PERTURBATIONS 

In this section we briefiy review the gauge-invariant ap- 
proach to odd-parity gravitational perturbations pq . As 
an application we derive a coordinate- invariant version of 
the RW equation ||l^ . We finally recall the arguments es- 
tablishing the stability of the Schwarzschild metric with 
respect to vacuum perturbations. 



A. Background expressions 

We are analyzing odd-parity perturbations of spher- 
ically symmetric background configurations. A spheri- 
cally symmetric spacetime (M, g) is a warped product of 
M = Af/S0(3) and S"^ with metric 



g = g 



R^g. 



(1) 



Here g is the standard metric on S^, and g and R denote 
the metric tensor and a real-valued function, respectively, 
defined on the two-dimensional pseudo-Riemannian orbit 
space M with coordinates x", say. Here and in the follow- 
ing lower-case Latin indices (a = 0, 1) refer to coordinates 
on (M,g), while capital Latin indices {A = 2, 3) refer to 
the coordinates i? and ip on {S'^,g). The dimensional 
reduction of the Einstein tensor yields 

Gab = -^ (2RAR + (dR , di?) - l) gat - | VaV^T? , 



GAB^l[2RAR~R^R]gAB, 



(2) 



where the off-diagonal components vanish, GAb — 0. The 
operators with a tilde and the inner product ( , ) refer to 
the two-dimensional pseudo-Riemannian metric g, and R 
denotes the Ricci scalar of g. 



B. Coordinate-invariant amplitudes 

Arbitrary perturbations of spherically symmetric back- 
ground fields can be expanded in terms of spherical tensor 
harmonics. For odd-parity perturbations the transverse 
spherical vector harmonics, Sa = {*dY)A form a basis 
of vector fields on 5^, while the harmonics V^^S'^} = 
^{'^aSb + ^bSa) are a basis of symmetric tensor fields 
on S^; see Appendix O for details. (Here * denotes the 
Hodge dual with respect to the metric g, and the Y^^ are 
the scalar spherical harmonics, where the angular num- 
bers £ and m are suppressed throughout, i.e., Y = Y^"^, 
Sa = S^.) The odd-parity perturbations of g^^ are, 
therefore, parametrized in terms of a scalar field k and a 
one- form h = hadx'^ , 

Sgab^O, SgAb ^ HSa, SgAB ^"^kS/iaSb}, (3) 

where k and ha depend on the coordinates x^ only. 



A vector field X — X^d^ generating an infinitesimal 
coordinate transformation with odd parity is determined 
by a function f{x^), where 



X° = 0, X^ = f S^ = 



i?2 



^AB 



Sb. 



(4) 



Under coordinate transformations induced by X the per- 
turbations of a tensor field transform with the Lie deriva- 
tive of the corresponding background quantity with re- 



spect to X: St 



^u 



St 



^y 



Cxt 



Xf-iiv 



Using CxgAb 



SAR^Vb{R-y) and CxOAB = 2/V{^5s}, 
perturbations transform according to 



hb hb ~ f f 



the metric 



(5) 



In a similar way one obtains the transformation laws for 
the perturbations of the Einstein tensor. Also using the 
background properties GAb = and 2G^ ~ G^S^ one 
finds 



SGAb -^ SGAb + G j^Sb R V 
SGab — > SGab + G D^ {aSb}! 



(6) 
(7) 



C. Coordinate-invariant Einstein tensor 



The computation of the coordinate-invariant compo- 
nents SG™^ is considerably simplified by the follow- 
ing observation: In the gauge where the scalar ampli- 
tude K vanishes, henceforth called the off-diagonal gauge 
(ODG), the perturbation h coincides with the coordinate- 
invariant perturbation H defined in Eq. (JIQ). (It is ob- 
vious from Eq. (H) that the ODG always exists and fixes 
the gauge function / uniquely.) Hence, for ^ > 1, the cor- 
rect invariant tensors are obtained by computing SG"^J 
in the ODG, and by substituting H for h in the resulting 
expressions. For ^ = 1 all perturbations are off-diagonal 
anyway, and one obtains the correct expressions in terms 
of the invariant quantity d{R~'^h). 

It is a straightforward task to compute SGg i, in the 
ODG. Using the formulas ( |A3|) , (A4) and (A5) derived 
in Appendix IaI Eqs. (pi) and (19|) yield the expressions 



5G 



•inv I ^ I \-?a 

Ab \ODG - ^ < V 



i?4V[fc {K}R-^) 



A, 
^hb 



SGIT \oDG = 0, SG'Xb Iodg- ^{aSb} ^'h, (11) 



where 



One may now use the transformation laws for k and 
hb to construct the following coordinate-invariant com- 
ponents: 



SGab 



'Ab 



SGAb — hb G j^Sb 



and, for £ ^ 1, 



<5G™^ = SGab ~ kG ^V^^S'b} , 



(8) 



(9) 



We recall that the scalar amplitude n defined in Eq. (g) 
is not present for £ = 1, since then V^a^b} vanishes. 
However, by virtue of Eq. (M), this also implies that 
SGab is already coordinate- invariant. (In fact, SGab 
vanishes identically for £ = 1, as will be shown below.) 
Hence, for £ — 1 one needs only the invariant components 
defined in Eqs. (|g), which do not involve the amplitude 

K. 

As the (5G*™ are invariant under coordinate transfor- 
mations generated by X, the expressions (||) and (Q) 
will only involve coordinate-invariant combinations of the 
one- form h and the scalar k. In fact, for ^ ^ 1, 5G™^ 
can be expressed in terms of the manifestly coordinate- 
invariant one-form iJ, defined by 



H 



«'H^) 



(10) 



This definition is again limited to £ 7^ 1. For ^ = 1, where 
K is absent, we will see that the remaining perturbation 
via the invariant two-form d{R^^h) only. 



h enters JG™" 



X ={£-!){£ + 2). 



Here we have used the background property 2Gg = 
G^jS^g and the fact that SG'^^ = SGab in the ODG. 
Since hb coincides with the invariant amplitude Hb in 
the ODG, we may replace hb by Hb in the above expres- 
sions, which makes them manifestly coordinate-invariant 
for £> 1. For £ = 1 the second term in the expression for 
SG™f^ vanishes, and hi, appears only via the coordinate- 
invariant expression W[b{ha]R^'^)- We therefore end up 
with the manifestly coordinate-invariant expressions 



0, 



'AB 



-d^HVt.S 



{A^B} 



fl2) 



and 



JG^Yda:" = ^ {dt [i?M (R-'H)] + Xh] , (13) 



which are valid for all values of £, provided that H is 
defined according to Eq. (|l^) for ^ > 1, and according 
to H = h for £ = 1. Here d^ = ?d? denotes the co- 
differential operator for p-forms on {M,g), e.g., d^iJ = 

-V'^Ha, {d^dH)b^2V''V[bHa]- 

The linearized Bianchi identity implies that the Ein- 
stein equation for SG™^ is a consequence of the equa- 
tion for 5G'']l^ . In fact, the first equation is the inte- 
grability condition for the second one, as is obvious for 
vacuum perturbations: Applying the co-differential to 
R^SG^Xb = yields d^iJ ^ 0, that is, SG^^ = 0. (For 
£ = 1 this integrability condition is void, in agreement 
with the fact that SG'^g vanishes identically.) 



D. Local uniqueness and linear stability of the 
Schwarzschild metric 

As an application we consider vacuum perturbations of 
the Schwarzschild metric. The relevant equation for the 
odd-parity sector was first derived by Regge and Wheeler 
|l7[ , and brou ght in a gauge- invariant form by Gerlach 
and Sengupta |l6| . A gauge- invariant approach which is 
based on the Hamiltonian formalism was given by Mon- 
crief |8|. 

The linear stability of the Schwarzschild metric follows 
from the dynamical behavior of vacuum fluctuations. In 
order to establish the local uniqueness property one also 
has to exclude all stationary perturbations other than 
the Kerr mode. While the stationary perturbations do 
not need to be normalizable, they are, however, subject 
to certain boundary conditions following from asymp- 
totic flatness and regularity requirements. Both station- 
ary and dynamical perturbations must be analyzed sep- 
arately in the sectors i > I and i = 1. 

The vacuum perturbations with odd parity are ob- 
tained from Eq. (Q), which yields 



1 



i?M 



H 

i?2 



-hA 



H 

i?2 



0. 



(14) 



This equation holds for all values of £ and comprises 
the complete information. The usual way to derive 
the RW equation from Eq. (|lj) is to decompose the 
one-from H with respect to Schwarzschild coordinates, 
H = Htdt + Hrdr, and to use the integrability condi- 
tion to eliminate Ht- This yields an equation for Hr 
alone, which is then cast into a wave equation for the 
function (1 — 2M/r)Hr/r. This can also be achieved in a 
coordinate-invariant way as follows: Using the integrabil- 
ity condition d' iJ = to introduce the scalar potential <& 
according to H ^ ?d(i?$), one may integrate Eq. (14). 



This yields Eq. (15) below for the potential <i> instead of 

Here we proceed in a different way, which is also 
coordinate-invariant. The basic observation is that in 
two dimensions the field strength two-form assigned to 
a one-form is equivalent to a scalar field. We therefore 
introduce the scalar field ^ according to 



* = R^id ( 4t 



i?2 



where the factor R^ turns out to be convenient. Apply- 
ing the operator ?d on Eq. ( [l4| ) and using the above 
definition yields the wave equation 



-^ + ^^15) + !^ 



* = 0, 



(15) 



where the two-dimensional Laplacian of a function is 
A^ = — dM^, and where we have used *dd' = dM?. 
Equation (|lj) is the coordinate-invariant version of the 



RW equation. In fact, it generalizes the RW equation, 
since it is not restricted to perturbations of static back- 
ground configurations. (The fact that the RW function 
^ = R^*d{R~^H) and the scalar potential <&, defined by 
H — ?d(i?<I>), satisf y the same equation will be explained 
at the end of Sect. 1vb| .) 

The positivity of the RW potential for i y^ 1 follows 
from the general expression (g) for Gab, which yields 
the coordinate-independent vacuum background equa- 
tion RAR + {dR, dR) = 1. By virtue of this, Eq. (|l|) 
assumes the form 



'A^A + Vbw * = 



with 



V, 



RW 



N 
R^ 



[3{N - 1) + i{£ + I)] 



and N = {dR, dR). Hence, Vhw is positive for finite 
values of R if dR is space-like and £>2. 

We may now use standard Schwarzschild coordinates 
r and t, defined by 



Rir, t) = r. 



g = -A^5W + ldr^ 



(16) 



to cast the RW equation into its well-known form. For a 
Schwarzschild background with mass M we have N{r) = 
1 - 2M/r, S{r) = 1, A^ A = -d^ + NdrNdr, and thus 



dt^ 



drl 



S(^(^+l) 



"-T) 



* = 0, (17) 



with dr* = A^~^dr. For £ > 2 the potential is non- 
negative in the domain of outer communications, and 
vanishes only asymptotically. Therefore, Eq. ( [l7| ) ad- 
mits no unstable dynamical modes. Furthermore, well- 
behaved stationary modes with £ > 2 can also be ex- 
cluded in a rigorous manner by applying the argument 
given in Appendix H. 

It remains to discuss the perturbations with £ = 1, for 
which Eq. dl3) is immediately seen to admit the solution 
\/ R. Since A = 0, we may also directly integrate Eq. 
(p^), which yields 

(H\ 6M 

where 6aM is a constant of integration. At this point it 
is important to recall that for ^ = 1 the one-form H = h 
is not coordinate-invariant, but transforms according to 
H -^ H + R^d{f/R'^). This imphes that the solution 
of the homogeneous part of the above equation is a pure 
gauge. Hence, with respect to Schwarzschild coordinates, 
the only admissible solution of the perturbation equa- 
tions (stationary and non-stationary) is H = 2a{M/r)dt. 
Using 5^"^ — and S^^ — — sin^i^, one finds with Eq. 
(I) 



ogtio = -a sm if , 

r 

which describes the Kerr metric in first order of the ro- 
tation parameter a. In conclusion, we have established 
the well-known result that the only physically admissi- 
ble odd-parity vacuum perturbation of the Schwarzschild 
metric lies in the sector £ = 1 and describes the station- 
ary Kerr mode. 



R "^{w^ — l)Tr?l, using the fact that dr^ commutes with 
*drr, and [dTr,Tr] = ^dr^. One finds 



Aw — 



1 



i?2 



(21) 



where d^ = ?dS, and Aw = — dMw — V^VaW. The 
Einstein equations, Gnu = STrGT^^, are obtained from 
the formulas ^ and (po|). Also using T^ = 0, one finds 



III. PERTURBATIONS OF YANG-MILLS FIELDS 

We are interested in perturbations of spherically sym- 
metric EYM solitons and black holes which give rise to 
odd-parity metric excitations. Before deriving the gauge- 
and coordinate- invariant expressions for the stress-energy 
tensor and the YM equations, we briefly recall some fea- 
tures of the background configurations. 



A. Einstein- Yand-Mills Background configurations 

The spherically symmetric EYM background configu- 
rations are assumed to be purely magnetic fl^, but not 
necessarily static. The metric is given by Eq. (||), while 
the gauge potential is parametrized in terms of a scalar 
field wix'') on M, 



^ = (1 — w)SdTr 



(18) 



where r^ = r • e^. Here the Tk = CTfc/(2j) are the su(2) 
generators, e^ is the radial unit vector in IR , and the Ufe 
are the constant Cartesian Pauli-matrices. The total ex- 
terior derivative of the vector valued function e^ is 6''^e^ 
(with A — -d^Lp), implying that 

dr^ — T§d{) + Tip siniJd^?. 

(See Appendix O for details.) Since r^ is an eigenfunction 
of the spherical Laplacian didr^ = — 2Trdri, the back- 
ground field strength, F = dA + A A A, becomes 



F = -dw A idTr + {w^ - l)TrdQ . 



(19) 



Using this expression, the components of the stress- 
energy tensor, T^, = .LTr {F^„F," - ig^.F^^F"^}, 
with respect to the background metric IM become 



Tab 

Tab 



1 



47ri?2 

1 
47ri?2 



1~ ^o . (w2-l)2 
2WaWh - -gab I 2WcW H — 



9AB- 



j w^ - 1)2 
2i?2 



TAb = 0, 



(20) 



where Wa = VaW, and where Tr { } denotes the normal- 
ized trace, Tr {^f} = 1. 

The background YM equation, D*F = d*F+[A, *F] = 
0, is obtained from the expression *F = —idw A dr^ + 



29ahR 



l^a^bR 



r ^ 



2WaWb 



l9ab y 



2WrW^ 



fl2 



l-iA(i?2) = G 



— W~ 



(22) 
(23) 



Equations ( |2l| ) - (|2^) are the spherically symmetric EYM 
equations in coordinate-invariant form. In the static case 
we may evaluate these expressions for the metric dlq), 
which yields (a prime denoting the derivative with re- 
spect to r) 



\ {NSw')' = w^ 



1 



(24) 



for the YM equation (pl|), and, with N{r) = 1 — 2m(r)/r, 



G 



(w2 - 1)2 



m 
b r 



l\2 



2N{w') 



(25) 
(26) 



for Eq. (|23| ) and for the trace- free part of Eq. (|22| ), 
respectively. Two special Abelian solutions to Eqs. ( |24| ) - 
( Pq ) are the Schwarzschild metric, m{r) — M = constant, 
5 = 1,^ = 1, and the RN metric with mass M and unit 
magnetic charge, A^ = 1 — 2M/r + G/r"^, S —\, w — 0. 

Asymptotically flat non- Abelian solutions with finite 
energy and nontrivial gauge fields are the solitons found 
by Bartnik and McKinnon R], and the corresponding 
black holes with hair g. They are obtained by numeri- 
cal methods and by analyzing the local solutions at the 
singular points of Eqs. (24)-(p6|), that is, at the origin, 
r = 0, the horizon, N{rH) — 0, and at infinity, r = oo. 
The local background solutions are given in Appendix 
y, since their behavior will be crucial to the existence of 
regular singular points of the perturbation equations. 



B. Gauge- and coordinate-invariant Yang-Mills 
perturbations 

In Appendix O we construct a convenient basis of 
su(2)-valued spherical harmonic one-forms. The odd- 
parity perturbations of the YM potential are then given 
in terms of two one-forms, a and /?, and three scalar 
fields, ^, v and a, over M, 

^^(<'>i) ^ Xia + X2f3 + fiTrdY + vYdTr + o-VATa, (27) 



where Xi, X2, and X^ are a scalar basis of su(2)-valued 
spherical harmonics, 

X^^YTr, X^^g^'^TAVBY, X3 = v^''taVbY, 

while Y = y^™ denote the ordinary spherical harmon- 
ics. (The antisymmetric tensor fjAB is defined by *d = 



fjie^ 



As usual, the cases i 



1 and £ — must 
be treated separately: For £ — 1, one has VX2 = 

—Y^^^^'d-Tr, implying that v and a combine to a single 
amplitude. Hence, 



(5^(fci) ^ Xia + X2P + flTrdY + vYdTr 



(28) 



In contrast to the gravitational sector, the odd-parity 
YM sector is not empty for £ — Q. As y^^^'') is constant, 
5 A is parametrized in terms of the one-form a and the 
fmiction i/, 



M(^=o) 



vdTr 



(29) 



One may now study the behavior of 5 A under gauge 
transformations, 5 A -^ 5 A + Dx, and under coordinate 
transformations, 5A -^ SA + CxA. Here D is the gauge 
covariant derivative with respect to the background con- 
nection (18), X is an su(2)-valued scalar field with odd 
parity, and Cx is the Lie derivative with respect to the 
infinitesimal vector field X defined in Eq. (Q). Con- 
sidering both gauge and coordinate transformations, the 
following results are established in Appendix pt 

For £ > 1 the metric perturbations are originally 
parametrized in terms of the function k and the one-form 
h, while the YM amplitudes are given by two one-forms, 
a and f3, and three functions, /i, v and a. Using the 
complete gauge and coordinate freedom, the entire set 
of perturbations reduces to three one-forms, H, A and 
B, and one function, C, all of which are invariant under 
both coordinate and gauge transformations. Adopting 
the ODG (k = 0) and the YM gauge fi = a = 0, the 
quantities H, A, B, and C, coincide with the original 
amplitudes h, a, (3, and 1^. [See Eqs. (|e|) and ( E12 ).] 
Hence, all physically relevant perturbations with £ > 1 
are given by 



a and b, defined in Eq. 
therefore given by 



(E4). The perturbations are 



Jab — 



Sg': 



{1=1) 



= 0, 



JAb — 



hbSA, 



M(^=^) = Xia + X26, 



(31) 



where a and b are gauge- invariant, but neither the metric 
nor the YM perturbations are invariant under coordinate 
transformations. The linearized EYM equations involve, 
however, only the gauge- and coordinate-invariant com- 
binations 






b + y^J^, 



(32) 



and d(R^^h), as we shall see later. 

For £ ~ there exist no metric perturbations in the 
odd-parity sector, and the YM perturbations are com- 
prised within a single gauge- invariant one-form a, defined 

in Eq. (Esl), 



6gl!='^=0, M(^=o)^ 



(33) 



IV. THE PERTURBATION EQUATIONS 

In this section we give the equations governing the odd- 
parity perturbations of a spherically symmetric soliton 
or black hole EYM background configuration. The am- 
plitudes are parametrized in terms of the gravitational 
one-form H, the YM one-forms A, B, and the YM scalar 
C. All amplitudes are gauge- and coordinate-invariant 
and, as we are not introducing specific coordinates, the 
resulting equations are not limited to static background 
configurations. The derivations are considerably simpli- 
fied by adopting the ODSG and by taking advantage of 
the su(2) harmonics constructed in Appendix |D|. How- 
ever, as the computations are still length y, w e disc uss 
only the basic steps in sections [V A , IV B and |IVC for 
i > I, i = 1 and £ = 0, respectively, and give a self- 
contained compilation of the results in Sect. [VD. 



^9ab 

j^(^>i) =XiA + X2B 



e^' - 0, 



^9Ab ^ HbSA, 

CYdTr, 



(30) 



with gauge- and coordinate-invariant amplitudes H, A, 
B, and C. The ODG for the metric perturbations, to- 
gether with the YM gauge fi = a = will be called the 
off-diagonal standard gauge (ODSG) henceforth. In the 
ODSG all gravitational and YM perturbations coincide 
with the corresponding coordinate- and gauge-invariant 
quantities. 

For £ = 1 the metric perturbations are already off- 
diagonal and there exists a gauge for which the YM 
scalars /x and v vanish, and the remaining amplitudes, a 
and /?, coincide with the two gauge-invariant one-forms 



A. Equations for £ > I 

For ^ > 1 we may proceed in the ODSG for which 
the metric and the YM perturbations coincide with the 
gauge- and coordinate- invariant amplitudes H, A, B, and 
C: 



^^(^>l) ^ XiA + X2B + C YdTr 



(34) 



We start by computing the coordinate-invariant stress- 
energy tensor. According to Eqs. (g) and (g) we have 



O^ab —°^ab i 



-AB 



Crr^ODG 



(35) 



and 






(36) 



since k = and Ha = ha in the ODG. The ST^i, consist 
of perturbations arising from variations with respect to 
the metric and the YM fields, 5Tf^i, 
where 



OgJ^fj.!/ + OA-Lfj,!^, 



^,T,. = -^Tr{iF„,^" 



/3> 



Jtii^ 






1 



TffA'i/ 



7 



5gap^ 



and 



In the ODSG the hnearized field strength, 5F = USA, 
is obtained from the formula (34) for 6A'^^^^\ Recalling 
that D is the gauge covariant derivative with respect to 
the background potential (|lj), one finds, also using the 
identities (D5), 

j^(^>i) = XiAA + XadS -X3,CAVL-B h VX2 

+ {wB - A)ATrdY + {dC~wA) AYdTr. (37) 

Using this, as well as the expression (nm for the back- 
ground field strength F and the formulas ( |34| ) for the 
metric perturbations, we end up with 



sz 



mv 
b 



0, ST 



mv 
AB 



1 

47r 



[B,&w)2V^aSb}. (38) 



and 



"-l^Ab 

Sa 
47ri?2 

{H , Aw)wi, 



(«;2 



l){Ai,~wB,,) 
{dw , dw) + 



+ R^{dB)baW''^Cwb 

(«;2 - 1)2 



i?2 



Hh 



(39) 



where we recall that all amplitudes are gauge- and 
coordinate-invariant. Here and in the following we use 
the obvious notations w"' = g°'^VbW and ( , ) for the inner 
product with respect to g, e.g., {H , dw) = g°'''Ha'^bW. 
[There is no factor 1/2 in front of the last term in Eq. 
(p9|), since, according to Eqs. (EOt) and (pq), ST^f^^ and 
6TAb differ by the term {Sny^R^iw^ -lyHbSA in the 
ODG.] 

The Einstein equations, SG'^J' = SttGJT;™, are now 
obtained from the above expressions and the formulae 
(|l|) and (|l|) for SG'^ . We have already argued that 
the (Ai3)-equation, 



d^H = -AG (B , di 



(40) 



is a consequence of the (y46)-equations and the linearized 
Bianchi identity. While this was obvious for vacuum per- 
turbations, one now needs the YM equations given below 



to verify this fact. Hence, the only independent Einstein 
equation is the one for the coordinate-invariant metric 
one- form H , 



jt 



RMt§ 



\H 



= 4G(w2 - 1) 



A-wB- 



w'-l 



i?2 



-H 



AG [iiR^dB + dwA H)idw - Cdw] , 



(41) 



where we also recall that A = {£—l)(i + 2). Here we have 
used the identities {H , dw)dw — {dw , dw)H — i{dw A 
H)idw and {dB)abw''dx°' = {idB)idw. 

The linearized YM equations also involve perturba- 
tions of both YM and metric fields. The latter arise from 
the variation of the Hodge dual in S(D*F) = 0, and yield 
the terms on the RHS of the following general expression: 

D*6F+ [SA, *i^] = 2 D * J^ - d ( -^(=3- I A *F, 



where JF^i, = F,"5gjj^„. Since the dual of this is an equa- 
tion between one-forms, and since the odd-parity basis 
of one- forms is five-dimensional for £ > 1, we obtain five 
equations. Again, the computation is considerably sim- 
plified in the ODSG for which we may use the gauge- 
invariant perturbations given in Eqs. (p3). As expected, 
it turns out that two YM equations can be obtained from 
the remaining ones. Using the tools developed in Ap- 
pendix |D|, we eventually end up with the following set of 
equations for the one-forms A, B and the scalar C: 



dt {R^dA) + [A + 2(w2 + l)]A - 2[A + 2]wB 

u;2-l 



2wdC + 2Cdw = (A + 2)- 



i?2 



-H. 



(42) 



d^ {R^dB) - 2wA + [A + (w^ + i)]B + dC 

^2- 1 
d\HAdw)-w „„ g, (43) 



i?2 



C = R^d^B- {dw, H) 



(44) 



The remaining two YM equations are the integrability 
conditions for Eqs. (|2|) and (|3|). Also using Eq. (|4|), 
these become 



dt 



A- 



l-w-" 
~R2- 



-H 



= -2{B,dw), 



(45) 



and 



AC-[A + (w;2 + l)] 



C_ 

i?2 



= 2(A, dw) -wd^A+[X + 2\ 



{dw , H) 

^2 



(46) 



Since Eqs. (^), ( [451 ) and (|4^) are consequences of 
the remaining equations, the complete system of pertur- 
bation equations consists of the three coupled equations 
(pl|), (^) and ( p3|) for the three gauge- and coordinate- 
invariant one-forms A, B and H, where C is given by Eq. 
( [44| ) . It will also turn out to be convenient to write these 
equations in terms of the one- forms A and B, defined by 



A^A+§, B^B + u.§, 



(47) 



in terms of which Eqs. (41), (W2) and ( [43| ) assume the 
form 

dt {R^Fh) + XH 
= AG [{R'^iFB)idw - Cdw + {w^ -1){A~ wB)] , (48) 



d^R^FA) +\[A~2wB 



-2{w^ + 1)A + AwB - 2Cdw + 2wdC, 



dt (R^Fb) +\(b~ w^ 
2wA - (w^ + 1)B - dC, 



with 



C^R-" 



d^B ~ wd^ 



H 
R^ 



(49) 



(50) 



(51) 



Here we have introduced the two- forms Fa, Fb and Fh, 
which are defined in terms of H, A and B as follows: 



FA = dA-FH, FB=dB-wFH, 
H 



^.-d(^ 



(52) 



i.e., Fa = dA^FB = dB + dw A Rr'^H. The three equa- 



tions (Uq) -(50) for the invariant one- forms A, B and H, 
with C according to Eq. (^l|), govern all physical odd- 
parity perturbations with £ > 1. We shall now argue that 
these equations hold for ^ = 1 as well, provided that one 
sets C = 0. 



B. Equations for i — 1 

For i — 1 the metric perturbations are off-diagonal and 
described by the one-form h, while the YM potential is 
parametrized in terms of two one-forms a and b, 



^9m^^ = hbSA, <^A(^=i) = Xia + X^b. 



(53) 



Although a and b are gauge- invariant, they are not in- 
variant under coordinate transformations, and neither 
is h. As the linearized YM and Einstein equations are 
coordinate-invariant, these will only involve the gauge- 



and coordinate-invariant one forms a and b defined in 
Eq. (H). 

The perturbation equations for £ — 1 are obtained from 
the equations for £ > 1 as follows: The linearized field 
strength two-form, 6F ~ DSA, for the background po- 
tential (|lj) and the perturbation (|53| ) becomes 

5^(fci) = Xida + X2db 

+ {wb - a) A TrdY + {b - wa) A Ydrr . (54) 



Formally, this is also obtained from the expression (|37_ 
for SF^^^-^^f by substituting a for A, b for B and by setting 
C — 0, where one also has to use VXj ~ = — y^^^^Mrr; 
see Appendix O for details. Hence, the invariant stress 
energy tensor for ^ = 1 is obtained from the expressions 
( p8| ) and (J3|) for £ > 1 by applying these substitutions 
and by replacing h for H. This yields 



and 



JTXrdx^ 



Sa 



47ri?2 

Sa 
AttR^ 



[{w^ - l)(a - wb) + R'^{idb)idw] 

Sa. 



(w^ - 1)2 ■ 
?(dw A h)idw ^^t; h 



i?2 



The coordinate-invariance of the last expression becomes 
manifest by writing it in terms of the one forms a and b 
given in Eq. (p2[). One finds 



STZ'^dx" 



Sa 
47ri?2 

Sa~ 

— * 
An 



{w^ - l){d -wb) 



d6- 



^d(^ 



Jdw, 



(55) 



where the metric perturbation enters only via the invari- 
ant two- form d{R~^h). 

It is now obvious that the complete set of linearized 
EYM equations in terms of the gauge and coordinate- 
invariant amplitudes a and b is obtained from Eqs. ([48|)- 
(pO|), by substituting A ^ a, B -^ b and C -^ 0. As we 
also have to substitute H ^ h, the LHS of Eqs. ([48|)- 
( pO| ) would, at a first glance, involve the non-coordinate- 
invariant amplitude h. However, since A = for £ = 1, 
the terms involving h itself vanish identically. We also 
point out that the algebraic equation ( pl| ) for C is not 
present ioi £ = 1, because the basis of one- forms is re- 
duced by one dimension. The complete set of pertur- 
bation equations in the sector £ = 1 thus assumes the 
surprisingly simple form 

dt {R^Fh) -AG [{R'^*Fb)idw + {w^ - l)(a - wb)] = 0, 

d^ {R^Fa) +2{w^ + l)a - Awb = 0, 

d^ (R^Fb) -2wa + (w^ + 1)5 = 0, (56) 

where a and b are the gauge- and coordinate-invariant 
one- forms given in Eq. (p3), in terms of which the two- 
forms Fh, Fa and Ff, are defined by 



F^ = d ( — I , Fa = da~Fh, Fb = db- wFh. 



C. Equations for £ = Q 

For a — Q there exist no metric perturbations with 
odd parity, and 5 A can be expressed in terms of a single 
gauge-invariant one-form 

The field strength is obtained by setting Y — \, B — Q, 
C = in the expression (p7| ) for 5F'^^^^\ and by substi- 
tuting a for A, 

Sp{i=o) ^Trda-wahdTr. 

The correct perturbation equation is now obtained from 
Eq. (H) by setting B = C = H ^ 0, where A = -2 for 
£^0. Also substituting a for A, Eqs. (Egl) and (ph yield 



d^ (R^Fa) + 2'u?a = 0, with Fa = da. (57) 



D. Summary 

All odd-parity perturbations of spherically symmetric, 
not necessarily static EYM configurations are governed 
by the three equations 

dt {R'^Fh) + XH 
= AG [{R^iFB)iAw - CAw + (w^ -1){A- wB)] , (58) 



d^ (R^F, 



X{A~2wB + w^^^ 



-2{w^ + l)A + AwB - 2CAw + 2wdC, 



dt {R^Fe 



X\B-w 



H 

R^ 



2wA ~ (w^ + 1)B - dC, 



(59) 



(60) 



for the three gauge- and coordinate-invariant one-forms 
A, B, and H, where X={£- !){£ + 2), 



C = (1 - 5^ - Si)R^ 



d^B - wd'' 



H 

i?2 



(61) 



and 



H 

i?2 



FH = d{-^], FA = dA-FH, FB = dB-wFH. 



The above equations are valid for all values of £, where 
only Eq. (p9h with H = B = is present for £ = 0. 



However, the expressions for the gauge- and coordinate- 
invariant amplitudes in terms of the original metric and 
YM perturbations are different for £ > 1, £ = 1 and £ = 0, 
respectively; see Appendix ^ 

£ > 1: The original metric perturbations are described 
by a one-form h and a function k, while the YM pertur- 
bations are given in terms of two one-forms, a and f3, 
and three functions, /i, v and a: 

Sgab = 0, SgAbdx'' = HSa, Sqab ^2n\/ s^aSb}, 
6 A = Xia + X2P + ^iTrdY + vYdTr + a^Xi. 

In terms of these amplitudes the invariant quantities ap- 
pearing in Eqs. ([5^)-(60) are, according to Appendix ^, 



h ( 2 '* \ 

— -d\ii^wa^w—). 



A = 



S^/? + u.^-d^. + z.^J 



H 



\\{^\ 



R^d 



\R^J 



(62) 



£ = 1: The original metric perturbations are described 
by the one-form h, while the YM perturbations are given 
in terms of two one-forms, a and /?, and two functions, 
fi and V. 



Sgab = 0, SgAbdx^ = hSA, Sg 



AB 



0, 



5A = Xia + X2I3 + ^TrdY + vYdTr. 
The invariant quantities now are 



(63) 



^""+i?^ 



H = h. 



ji — wv 
1 — w^ 

W/i — V 



(64) 



i = 0: There exist no metric perturbations, and the 
YM perturbations are given in terms of a one-form, a, 
and a function, v. 

5 A = TrOl + U dTr. 

In terms of a and v the invariant quantity A is given by 
A = a-d{-), (65) 



and, as mentioned above, the correct perturbation equa- 
tion is Eq. (^ with A = -2, C = 0, B = and iJ = 0. 



V. NON-ABELIAN STABILITY AND LOCAL 

UNIQUENESS OF THE 

REISSNER-NORDSTROM SOLUTION 

For w = the static, spherically symmetric EYM equa- 
tions (E4h-(E6h admit the RN solution with unit magnetic 



charge. The stabihty and local uniqueness properties of 
the RN metric with respect to non-Abelian perturbations 
are, therefore, obtained from Eqs. (pq)-(|60|), which de- 
couple into two sets for w = 0: The first set, involving 
the one-forms H and A only, is obtained from Eqs. (pq) 
and (H), 



dt 
dt 






XH = -iGA, 



+ (A + 2)A = for e>l. (66) 



Since w = 0, the remaining equation for B does not con- 
tain the amplitudes H and A. Using C — R^d'B for 
£ > 1 and C = for £ = 1, we have 

d^ {R^dB) + d (R^d^B^ + {X + 1)B = for ^ > 1, (67) 
d^ (R^dB) + B ^ ioT e = 1. (68) 

Since A is the gauge- and coordinate-invariant version 
of the amplitude in front of the isospin harmonics TrY^, 
Eqs. (pq) govern the Abelian part of the perturbations, 
that is, Einstein-Maxwell perturbations of the RN metric. 
In contrast to this, Eqs. ( |67| ) and ( |68| ) for B are not 
present in the Abelian case, and describe non-Abelian 
perturbations of the RN metric with i > 1 and £ = 1, 
respectively. 

A. Perturbations with I > 1 

We start with Eqs. (£^) describing the Abelian part of 
the perturbations. For £ > I the integrability conditions 
for these equations are d'A = and d' H = 0, implying 
the existence of two scalar fields, '^h and ^a, defined by 

Sd (i?*H) = VXH, id^A = ViGA. 

Substituting ^/f and ^a for H and A in Eqs. (p9), and 
integrating both equations yields the following coupled 
wave equations for the scalar fields $// and ^^: 






dR 



Rd'r^ +-. 



i?2 



A 



yiGX 
~1^ 



* 



H 



i?2 

A + 2 



ViGA^ 



4G 



i?3 



*. 



For w = the background equation (|2^) becomes 
i?3dt(di?/i?2) = 3(di?, di?) - 1 -H G/R"^. Using this and 
introducing standard Schwarzschild coordinates, R = r, 
{dR,dR) ^N ^l- 2M/r + G/r'^, yields 






A + 2 + ^r 



4G\ 



J_ / -3M V4GA \ 
'r3 iv\/4GA 3Af j 






= 0. (69) 



The above equation was first obtained by Moncrief by 
different means |1^. Since the off-diagonal part of the 
potential is symmetric and constant, Eq. ( |69| ) can be 
decoupled. Using the non-negativity of N{r), as well as 
the regularity condition M > G, the eigenvalues of the 
potential are found to be positive, implying the absence 
of unstable modes. Taking advantage of the argument 
presented in Appendix |B| stationary modes are excluded 
as well. (The eigenvalues of the potential are positive for 
finite r and behave like £{£ — l)r~^ + 0{r~^) for r — > oo, 
implying that the asymptotically finite solutions behave 
like r~^ .) Hence, there exist neither unstable modes nor 
admissible stationary solutions to Eqs. (66) for ^ > 1. 



In order to discuss the non-Abelian part of the pertur- 
bations we introduce the scalar fields Hi = R^idB and 
112 = R^d^ B. In terms of these, Eq. ( [6^ ) assumes the 
form 



fdHi + dn2 + (A + l)B = 0, 



(70) 



which can also be viewed as the Hodge decomposition of 
the one- form B (see the comments below) . Applying the 
operators ?d and d = ?d? to this, it is immediately seen 
that Hi and H2 are subject to the same equation, namely 



AH,; 



A±ln 



1,2, 



(71) 



where we recall that A = {£~l){£+2). With respect to the 
static RN background, R{r,t) = r, g ^ N{~dt'^ + drl), 
with N{r) = 1- 2M/r + G/r"^ and dr^ = dr/N, one has 



dt^ 



drl 



N{r) 



£U+l) 



H, = 0. 



(72) 



Since the operator is positive, unstable modes are absent. 
Furthermore, well-behaved stationary solutions are ex- 
cluded as well, since the potential is of the type required 
to apply the argument given in Appendix H. [Also note 
that the one-form B is obtained directly from Hi and H2 
by Eq. (^.] 

As we shall continue to use the above method, it is 
worthwhile noticing the following: In two dimensions an 
arbitrary one-form 0, say, gives rise to two scalar fields, 
gi = d'6' and 52 = *d6. On the other hand, the Hodge 
decomposition of a one-form in two dimensions involves 
two different scalar fields, 9 = d/i + ?d/2 (provided that 
the harmonic part vanishes). If 9 is subject to a linear 
wave equation, then the latter gives rise to an algebraic 
relation between the two different parameterizations, al- 
though the two scalar pairs are defined on different dif- 
ferential levels. (This is also the reason why, in Sect. 



II D, we have obtained the same RW equation (15) for ^ 
and $, defined by ^ = RHd{H/R^) and H = Sd(i?$), 
respectively.) 
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B. Perturbations with I — 1 

Defining Hi = iRMB as for ^ > 1, Eq. (|8|) for B 
reduces to *dUi + B ^ 0. Applying the operator ?d 
yields the same equation for Hi as before, that is, Eq. 
( jn\ ) , where now i = 1. As the potential remains positive 
for ^ = 1, we conclude that Eq. ( |68|) admits neither 
unstable modes nor admissible stationary perturbations, 
which establishes the stability and the local uniqueness 
of the RN metric with respect to non-Abelian odd-parity 
perturbations. 

It remains to consider Eqs. ( |6q ) for £ ~ 1, i.e., for 
A = 0. As these equations are also present in the Abelian 
case, we will recover the absence of unstable modes, while 
the only stationary perturbations are those describing 
the Kerr-Newman excitations of the RN solution. This is 
seen as follows: For A = the only integrability condition 
for Eqs. ( p6| ) is d' A = 0. Using this to define the scalar 
field ^ according to 



Eqs. 





Sd* = A 


) 


(73) 


can be integrated, which yields 






(74) 


B? 




- 2* == 0, 


(75) 



where QMa is a constant of integration, and where we 
have used the fact that ^ is defined up to a constant 
in order to neglect the second constant of integration. 
Eliminating the gravitational perturbation H from the 
above equations yields the following inhomogencous wave 
equation for VP: 



-A- 



2 AG' 



* 



6Ma 



(76) 



As the operator on the LHS is positive, we conclude 
again that there are no unstable modes. Using standard 
Schwarzschild coordinates, R — r, N ^ l — 2M/r + G/r'^, 
we have A = —N~^df + drNdr, and the inhomogencous 
problem admits the particular solution '^ — a/r. By 
virtue of Eq. (|74| ) and definition ( |7^ ) this yields, up to a 
gauge. 



N 
H^a{N-l)dt, A = a—dt. 



(77) 



Recalling that for £ = 1 one has 5ga^=Q, Sgaip — 
-Ha sin^i? and SA= {A- H/P)Xi + {B ~ wH/r^)X2, 
we find with w = and B = 



Sg, 



tip 



SA: 



f2M G 
\ r r^ 

-Tr cosz9di, 



sin^??, 



which is the Kerr-Newman excitation of the magneti- 
cally charged RN metric. In order to see this, we com- 
pute the electric field, 5E = -5F{dt, ■) = -T)5A{dt, ■) = 
aTrd{cos-d / r"^) , where we have used DTr = and w = 0. 
Hence 



5E 



= —T^a I 



sin^ 



Tj-a 



dd^ 



2cos?9 



-dr 



Since the magnetic field of the background solution is 
B = ~Tr{*dn){dt,-) = -Tr(l/r2)dr [see Eq. (^ for 
w = 0] , we obtain indeed the magnetically charged Kerr- 
Newman solution in first order of the rotation parameter 



C. Perturbations with £ — 

Since the odd-parity gravitational sector is empty for 
i = 0, the perturbations of the RN solution are governed 
by Eq. (|9|) with H = B = 0, w = and \ = -2, 

di (R^dA) = 0. 

With respect to Schwarzschild coordinates the solution 
is A = {q/r)dt, where g is a constant of integration. The 
perturbation of the gauge potential now becomes SA — 
Tr{q/r)dt, which gives rise to a radial electric field. 



6E 



-Tr—dr. 



Hence, we obtain the embedded magnetic RN solution 
with infinitesimal electric charge q. (Note that the metric 
remains unchanged in first order of q.) 

In conclusion, we have shown that the RN solution 
is stable with respect to both Abelian and non-Abelian 
odd-parity perturbations for all values of ^. Also, the only 
physically admissible stationary modes are the Abelian 
ones, describing electric Kerr-Newman (^ = 1) and RN 
{£ = 0) excitations of the magnetic RN metric. 



VI. NON-ABELIAN STABILITY AND LOCAL 

UNIQUENESS OF THE SCHWARZSCHILD 

SOLUTION 

The Schwarzschild metric solves the spherically sym- 
metric EYM background equations with w = 1. As the 
stress-energy tensor is quadratic in the field strength, the 
gravitational perturbations decouple in first order for all 
values of £, and are governed by the RW equation for 
vacuum perturbations. The remaining equations, de- 
scribing Abelian and non-Abelian perturbations of the 
Schwarzschild metric, admit no unstable modes, and, for 
£ > 1, no acceptable stationary excitations either. For 
i = 1 the only stationary YM perturbation is the RN 
mode in the Abelian sector. 
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A. Perturbations with £ > 1 

The gauge- and coordinate-invariant one-forms A, B 
and H given in Eqs. (p2) for £ > 1 are well-defined for 
w ~ 1. The perturbations are, therefore, governed by 
Eqs. (pSl)-®, where Eq. dssl) decouples for w = 1 



and reduces to the usual equation describing the vacuum 
perturbations of the Schwarzschild metric. 



^t 



rM 



Ry 



XH = 0. 



(78) 



In Sect. II D we have already recalled that this equa- 



tion admits neither unstable nor well-behaved stationary 
solutions for £ > 1. 

In order to discuss the non-vacuum perturbations of 
the Schwarzschild metric, it is more convenient to re- 
sort to the original one-forms A = A — H/R^ and 
B = B — wH/R'^, used in Sect. |^ to derive the per- 
turbation equations. In terms of A and B, Eqs. ( |59| ) and 
( |60| ) become for w = 1 

d^ {R^dA) - 2d (R^d^B) + (A + A)A - 2(A + 2)B = 0, 

d^ (R^dB) + d (R^d^'B) - 2.4 + (A + 2)B = 0. (79) 

The above system is equivalent to four coupled equations 
for four scalar fields. In order to decouple these equations 
completely, we note the following: The terms with B and 
d'B can be eliminated, which shows that the intcgrabil- 
ity condition is d A = 0. Using this, and applying the 
co-differential operator on cither of the above equations, 
yields a wave equation for the scalar field d* B alone. 



A + 2 



n. 



0, Hi 



R^d^B. 



(80) 



Since the integrability condition implies that the scalar 
n^ = R^d^ A vanishes, it remains to find the equations 
for the field strengths dA and di?, or, equivalently, for 
the scalar fields ^a and '^b, defined by 



*A = Si?2dyl, *s = VX + 2iR'^dB. 

Applying the operator Sd on Eqs. 
system 



79) then yields the 



-A- 



1 
R^ 



A + 4 



-2VA + 2 

-2VAT2 A + 2 






0, 



(81) 



which can be diagonalized, since the potential is symmet- 
ric and constant. The eigenvalues are 



(^+l)(^ + 2) 
£{£-!) 



X + 3± V4A + 9 = 



Having solved Eqs. (|8^) and (|8l|), the expressions for 
the one-forms A and B in terms of the scalar fields are 
obtained from the original equations WW: 



A = --id ( *A . , 

A V x/AT2 



*? 



B = - 



1 



A(A + 2) 



;d 2*A 



A + 4 

VAT2 



*f 



Adn, 



Since the operators in Eqs. (pO) and ( |8l| ) are positive, 
we conclude, using the argument given in Appendix ^, 
that the Schwarzschild solution admits neither unstable 
nor stationary non-Abelian odd-parity modes with £ > 1. 



B. Perturbations with £ — 1 

For w — 1 Eq. (^) decouples for all values of £. The 
vacuum perturbations of the Schwarzschild metric with 
£ = 1 are, therefore, governed by Eq. ( [7q ) with A = 0. 
We have already recalled in Sect. [ID that this equation 



cannot give rise to unstable modes, while it admits the 
well-behaved stationary solution H — {2aM/r)dt, giving 
rise to the Kerr excitation of the Schwarzschild metric. 



A 2M . 2 

OQtui = —a sm V. 

r 



(82) 



In order to analyze the YM sector, we first note that 
the gauge invariant quantities introduced in Eqs. (E4) 
for £ = 1 are not well-defined if it; = 1. Hence, the £ = 1 
perturbations of the Schwarzschild background require a 
special treatment: For w = 1 and £ — 1 we define a, b 
and c in the same way as for ^ > 1, that is, by Eqs. (p^). 
Hence, a — a — d/i, b = j3 and c = v — 11, where the 
one-forms a, (3 and the scalars /i, 1/ parametrize 6A^^^^' 
according to Eq. (pq). Since the gauge fields vanish 
on the background, all YM amplitudes are coordinate- 
invariant, and it remains to consider their behavior under 



gauge transformations. By virtue of Eqs. ( |El| ) and ( |E2| ) 
c remains invariant, whereas a and b transform according 
to —^ a+d/2, b — + fe+d/2. Repeating the arguments given 
in Sect. [VB, the perturbation equations for w = I and 
£ = 1 eventually become 

d^ {R^da) - 2dc + 4(a - 6) = 0, 
d'' (R^db) +dc~2ia-b) = 0, 
R'^d^ (a - 6) + c = 0, 

where c, da, d6 and (a — b) are gauge-invariant. Sub- 
tracting the first from the second equation, and using 
the third one to eliminate c, we obtain an equation for 
the one- form {b — a). This is decoupled in the usual way, 
that is, by introducing two scalar fields according to 

ci = R'^d\b~a), C2 = R^id{b-a). 

Applying the operators d' and ?d on the equation for 
(6 — a) yields the following wave equations for ci and 02'- 



-A 



i?2 



ci = 0, 



-A + -^)c2 = 0. (83) 
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Since the operators are positive, we may use the stan- 
dard argument to conclude that Eqs. (|8^) admit nei- 
ther unstable nor well-behaved stationary modes. Hence, 
ci = C2 = 0, implying that a = b and c = 0. It 
therefore remains to solve d{R^*da) = for the gauge- 
invariant scalar field Sda. With respect to Schwarzschild 
coordinates, the result is a = b = {q/r)dt plus gauge 
terms, where g is a constant of integration. Now using 
a — a + dfi, P — b and h' = c + fi in Eq. (Eq) gives 
j^(<;=i) ^ ^Xi + bX2 + cYdTr + d[iiXi), and thus, with 
c = 0, a = 6 = {q/r)dt and Xi + X2 = T2, 

6A^T^-dt 
r 

plus a pure gauge term. Using 5F — T)5A — dSA for 
w = 1, this gives rise to the electric field 



5E- 



-T^—dr. 



(84) 



The solutions (|8^) and ( |84| ) describe the Kerr-Newman 
excitation of the Schwarzschild metric in first order of the 
rotation parameter a and the electric charge q. 

C. Perturbations with ^ = 

The relevant perturbation equation is Eq. (B9|) with 
B = H = {), w = l and A = -2. The amplitude A == 
A is gauge- invariant and, by virtue of Eq. (pa), well- 
defined. Equation @) becomes d\R^dA) -I- 2A = 0. 
Using the integrability condition d^ A = 0, the scalar field 
^ is defined according to '5 = R^idA, in terms of which 
Eq. (KB) becomes 



(85) 



^A + ^)vI,.0, 



which admits neither unstable nor acceptable stationary 
solutions. (Note that the RN excitations with £ = of 
the Schwarzschild metric lie in the even parity sector.) 



In the stationary case, the excitations of a spherically 
symmetric EYM background decouple into two Sturm- 
Liouville problems, governing the electric and the mag- 
netic perturbations, respectively. The particular case 
a. — 1 was analyzed in Refs. pffl and m^ by different 
means. There we have shown that the electric sector gives 
rise to a two parameter family of slowly rotating and / or 
electrically charged black hole excitations, and to a one- 
parameter family of slowly rotating, electrically charged 
solitons. In this section we generalize these results as fol- 
lows: We show that for all values of ^ > 1 the electric 
perturbations are governed by a three-channel Sturm- 
Liouville problem, while the magnetic sector is described 
by a single Sturm-Liouville equation for £ > 1 and is triv- 
ial for € = 1. A careful analysis then reveals that neither 
the electric nor the magnetic sector admit well-defined 
stationary soliton or black hole excitations if ^ > 1. This 
establishes the result that the only stationary odd-parity 
modes of the BK solitons and EYM black holes are the 



ones found in Ref. |11| for ^ = 1. 

It turns out to be convenient to parametrize the two- 
dimensional background metric g in terms of the radial 
coordinate p, defined such that g becomes conformally 
flat. 



NS'^de + j-dr'' = a {-de + dp2 



N 



(86) 



with a{p) = N{r)S'^{r) and dr = NSdp. [The coordinate 
p generalizes the coordinate r* used in the Schwarzschild 
or the RN case. We also note that idt = — Sdp and 
(T?(d< A dp) = — 1.] The invariant one- forms A, B and H 
are expanded with respect to t and p, e.g.. 



H = Hodt + Hidp. 



(87) 



Since we restrict ourselves to stationary perturbations 
the coefficients Hq, Hi, etc. are functions of p only. As 
we shall argue below, the equations involving the zero- 
components, henceforth called electric perturbations, de- 
couple from the equations for the one-components, hence- 
forth called magnetic perturbations. 



VII. STATIONARY PERTURBATIONS OF 
NON-ABELIAN SOLITONS AND BLACK HOLES 



A. The electric sector 



Having analyzed the complete set of non-Abclian odd- 
parity perturbations (stationary and dynamical) of the 
Schwarzschild and the RN solutions, we now turn to the 
general case, that is, to non-Abelian perturbations of 
static non-Abelian background configurations. The dis- 
cussion of the corresponding perturbation equations is 
a considerably more involved task, since the techniques 
used above cannot be applied if w is not constant. Our 
primary goal in this section is to classify all stationary 
odd-parity perturbations of both the BK solitons 0] and 
the static, spherically symmetric EYM black holes pi . 



The electric perturbation equations involve the ampli- 
tudes Hq, Aq and Bq only. [For Z = 1, we may take 
Hq = ho, Aq = qq and Bq = bo, since, by virtue of Eqs. 
(g), (E4) and (Ell), these amplitudes are invariant under 
stationary coordinate transformations.] Using the fact 

that [d^(i?'^dA)]o = — 9p((T~^i?^(9pAo) for stationary per- 
turbations of a static background, the zero-components 
of Eqs. (|4l]), (III) and (H) may be cast into the following 
three-channel Sturm-Liouville equation: 



-dr^d + Kd- dK^ +L + P 



) II = 0, 



(88) 
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where v = (i/o/ViGr, Aq/C, -Bq j , C = VA + 2, and 
where the differential operator d is defined by 

f d _ f d 
a dp 5* dr ' 



with r and p according to Eq. (|8q). The 3x3 matrices 
K, L and i-* are given in terms of the background fields 
w, N and a = S'^N. The only non-vanishing matrix 
element of K is K13 — \fAGr dw, while the symmetric 
matrices L and P are 



(7 



and 



2N + X sym. sym. 

A + 2(f + u;2) syjn. | , (89) 

-2(w A + (l + w2) 




2Ga{dw) 



sym. sym. 
sy?7i. 




(90) 



The formally self-adjoint equation (pq) holds for all val- 
ues of £ > 1. [In particular, ioi £ — 1 it is equivalent to 
the Sturm-Liouville equation derived in ||20| , which was 
shown to admit the stationary modes mentioned above 
|ll| . However, the transformation between the two i = 1 
sets of equations is not algebraic, because the original 
formulation given in |20| was based on the generalized 
twist potential. 

Since Eq. (g8) has regular singular points at the origin, 
r = 0, at the horizon, r = ru (where N^rn) — 0), and 
at infinity, r = 00, it is possible to compute the number 
of stationary modes. Applying the standard theory (see, 
e.g., [gl[) we will now discuss the local solution spaces. 



2. The solution space at infinity 



The asymptotic expansions (C2) of the background 
quantities show that the leading order behavior of the so- 
lutions to Eq. (|88| ) is again completely determined by L: 
The solutions behave like r", where again a = — (^ + 2), 
-{£ + 1), -e, e - 1, e OT i + l. The space of asymp- 
totically flat solutions is, therefore, i/iree-dimensional for 
£ > 1, and /oMr-dimensional for £ = 1. For i = 1 the 
asymptotic expansion is found to be 



v{r) 



( ^^^ 
V r 

l + (l-7) 



O 



logr 



C3 



^) 



v% 


eo 




■)] 



-'^ 



(92) 



The constant C2 is proportional to the total angu- 
lar momentum (5J, while co and ci are proportional 
to the asymptotic value of the electric YM potential 
S^oQ and the electric YM charge 5Qe^ respectively: Us- 
ing the above expansion in the expressions (pi|) for the 
linearized local Komar integrals, we find [with *Fb — 
a-^B'o + w'Ho/R^) etc.]. 



SQeir 



00 



TCi, SJ(r 



00 



iOC2- 



Furthermore, the above expansion, together with the def- 
inition ( p3| ) and (5$ — SA{dt), shows that S^oo is propor- 
tional to cq. It is worthwhile recalling that, in contrast 
to the Abelian case, cq cannot be "gauged away". This 
is also obvious form the fact that the expression for SF 
involves an asymptotically vanishing term proportional 
to co/r, unless for w = 1. 



3. The solution space at the horizon 



1. The solution space at the origin 

The leading order behavior of the solutions to Eq. ( pq ) 
in the vicinity of the origin is determined by the cen- 
trifugal barrier L, as can be seen from the expansions 
(CI) of the background quantities. The solutions behave 
like r", where a = -{£ + 2), -{£ + 1), -£, £ - 1, ^ or 
£ + 1. Hence, the space of regular solutions at r = is 
t/iree-dimensional for all values oi £ > 1. The expansion 
becomes 



v{r) = dif 
+ d2r 



^+^^{2{£+2)b+l)br'e_+Oir') 
2b 



■0^2^+1 
+ d,/+'[e^+Oir)] 



re_+0{r^) 



(91) 



where Co = (1,0,0), e+ = (0,C, -0 and e_ = (0, C, ^ + 
1), and di, d2 and d^ are constants, and where b is the 
fixed constant appearing in the expansions ( plD of the 
background solutions. 



Using the background expansions (C3) at the horizon, 
the solutions to Eq. ( p8|) behave like (r — r//)", where the 
eigenvalues are a = and a = 1, and the multiplicity is 
three in both cases. For a = the three eigenvectors may 
pick up logarithmic terms in next-to-leading order, which 
destroy the regularity of the horizon. A careful analysis 
shows that the number of eigenvectors with logarithmic 
terms in next-to-leading order is equal to the rank of the 
symmetric matrix 



/A + 4G%:i 



Si = 




sym. 



\ 



sym. 
X + 2{l + wj^) sym. 



which is proportional to the leading order term of i + i-* 
in r — r/f . The determinant of Si is given by 

det 51 = A [A^ + (3 - wjjjX + 2{l-w%f + 8G Gj^] , 

where we recall that wh = wirn) and Gh = WHiwfj — 
l)/rH- This shows that the rank of Si is three ior £ > 1, 
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while one may also verify that the rank is two for £ = 1. 
Hence, all solutions with a = must be excluded, unless 
£ = 1, in which case there exists one acceptable eigenvec- 
tor. The physical space of solutions at r = r// is, there- 
fore, </iree-dimensional for £ > 1 and /owr-dimensional for 
£^1. 



4- Soliton excitations 

Since the BK background is continuous, and since the 
perturbation equations are linear with continuous coeffi- 
cients for < r < oo, the local solutions (^) and JQ^ ) 
admit extensions to the semi-open intervals [0,cxd) and 
(0, oo], respectively. Since, for £ — 1, these solution sub- 
spaces are three- and four-dimensional, respectively, and 
since the total space of solutions is sia;-dimensional, we 
conclude that the intersection space is generically one- 
dimensional. Hence, there exists (at least) one global 
solution, describing the rotating charged solitons found 

inQ. 

For £ > 1. the intersection space is generically trivial, 
since the solution spaces are three-dimensional at both 
the origin and infinity. Hence, there exist no generic soli- 
ton excitations for £ > 1. In fact, non-generic solutions 
are excluded as well, as we shall prove below. 



ri=diag(r,l,l), Ta 



[Note that the components of u 
amplitudes introduced in Eq. 





1/C 
w 

T^^v_ coincide with the 

H): VaGui = Ho/r^, 

= Bo + wHo/r^.] The Sturm- 

now assumes the desired form 

2rrT 



C,U2 = ^0 + Hq/t^ Us 

Liouville equation ( p8| ) 

(P3|), with the symmetric matrices A = r^T^ o T, 

S = tI 0S0T2, where 



S= - 

a 



Ar2 + 4G(w2-l)2 

\/4GC(l-w^) 
\/4Gw(u;2 - 1) 



sym. 

X+2{l + w^) 

-2Cw 



sym. 

sym. 

X+l + w'^ 



It is not hard to see that the matrix S is positive for 
all values of € > 1 and positive semidefinite for ^ = 1. 
Furthermore, by virtue of the expansions given above for 
£ > 1, the boundary term u ■ Adu vanishes at the origin, 
at the horizon, and at infinity. Both soliton and black 
hole solutions are, therefore, excluded as a consequence 
of the argument given in Appendix H. 

We emphasize that the boundary terms at the origin 
and at the horizon do give non- vanishing contributions 
if ^ = 1. The positive operator in Eq. (93) is, therefore, 
self- adjoint only for £ > 1. 



5. Black hole excitations 

Applying the same argument as in the soliton case, 
we conclude that Eq. (jSg) admits a two-dimensional in- 
tersection space of global solutions for £ — 1^ since the 
local solution spaces at the horizon and at infinity are 
four-dimensional. The solutions give rise to the black 
hole excitations found in [|ll| , which are parametrized by 
their total angular momentum 5 J an their electric YM 
charge 5Qe. 

For £ > 1 there exist again no generic solutions, since 
the solution spaces at the horizon and at infinity are 
three-dimensional only. It therefore remains to exclude 
non-generic solutions, which we shall do next. 



6. Absence of non-generic solutions for £ > 1 



Our aim is to show that Eq. mw with the bound- 
ary conditions discussed above admits neither soliton nor 
black hole solutions for £ > 1. We do so by casting Eq. 
( p8| ) into the form required to apply the argument out- 
lined in Appendix H. This is achieved by performing the 
linear transformation v = Tu, which yields 



{-dAd + S)u = {), 



(93) 



where A is symmetric and positive, while S is symmetric 
and positive semidefinite. The linear transformation T 
is given by T = Ti o T2 , where 



7. Conclusion 

We have proven the following local uniqueness theo- 
rems for odd-parity perturbations in the electric sector: 
The only stationary, asymptotically flat black hole solu- 
tions which are inflnitesimally close to the static, spheri- 
cally symmetric EYM black holes are the rotating and/or 
electrically charged excitations in the £ = 1 sector. The 
only soliton solutions which are infinitesimally close to 
the BK solitons are the electrically charged excitations 
in the £ = 1 sector. 

These results are in agreement with the non-Abelian 
staticity theorem |22], which asserts that spacetime is 
static and purely magnetic if the combination U,hJ — 
Tr{$tx)<5e} vanishes, where ^h is the angular velocity 
of the horizon: For / = 1, non-static solitons and black 
holes can exist, while, for £ > 1, there is no contribution 
to J and Qe [see the general formulae ( Fl )] , implying that 
the non-static and electric contributions Hq, Aq and Bq 
must vanish. 



B. The magnetic sector 

For stationary perturbations one has dr A = —a^^A[, 
where here and in the following a prime denotes differ- 
entiation with respect to the radial coordinate p, defined 
in Eq. (Pq). Since the background is static, one also 
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has {A, dw) — —a~^w'Ai. Hence, the gravitational con- 
strai nt (4C| ) and the YM constraints (^5|), ([46|), as well as 
Eq. (^) involve only the one-components of A, B and 
H. It is, therefore, possible to express the YM amplitudes 
Ai , Bi and C in terms of the gravitational perturbation 
Hi: 



Ai 



2. Black hole excitations 



Bi = 



C = 



(w^ 


_^^R^\Hi 
2g) i?2 ' 


AGw' '' 


1 
a 


,„ i?2 (H[\' 
AG \w' J 



(94) 



Using the above expressions and the circumstance that 
[(T {R^d{H / R'^))]i vanishes for stationary perturbations 
of a static background, the one-component of the gravita- 
tional equation ( [ll|) yields the following Sturm-Liouville 
equation for Hi: 

_d_J__d_ [£(i+l)-2w^]a-'iGi 
'd^^d^^ W 



w ■ 



Hi=0, (95) 



where we recall that w' = dw/dp = NSdw/dr. 

The above equation holds for £ > 1 only. For i = I the 
perturbations are governed by Eqs. (|5^). Since the one- 
components of the first terms in these equations vanish 
for stationary perturbations, we obtain ai = bi = 0, pro- 
vided that w'^ — 1 does not vanish everywhere. Now using 
the fact that there exists a gauge for which hp vanishes 
ii i = 1, we conclude that magnetic excitations cannot 
exist for £ = 1. (T he ca se w"^ — 1, £ — I has already been 
discussed in Sect. 



VIB.) 



Equation ( |95| ) has regular singular points at the origin, 
i? = 0, at the horizon, N = {dR , di?) = 0, at infinity, 
R = oo, and at all points where w' vanishes. (For the 
one-node background solutions this is only the case at 
the origin and at infinity.) In order to conclude that 
Eq. (p3) generically admits neither acceptable soliton 
nor black hole excitations, it is sufficient to discuss the 
regular singular points at the boundaries in leading order. 



1. Soliton excitations 



Using the expansions (CI) for the BK background at 
the origin shows that the fundamental solutions to Eq. 
( P5| ) behave like r^~^^ and r^~^. Since £ > 1, the subspace 
of solutions giving rise to finite metric perturbations is, 
therefore, one-dimensional at the o rigin . In the asymp- 
totic region one uses the expansions (C2) to conclude that 



the fundamental solutions behave like 



and 



implying that the subspace of bounded solutions is again 
one-dimensional. Generic soliton excitations are, there- 
fore, excluded. [The subspace of bounded solutions at 
the inner points w' — turn out to be two-dimensional. 
It is, however, generically not possible to match the so- 
lutions from r = and r = oo at the points w' = such 
that the amplitude G is continuous.] 



Using the horizon expansions (Cq) shows that the fun- 
damental solutions to Eq. ( p5| ) behave like (r — r//)° and 
(r — r//)2. The first solution is physically unacceptable, 
since the invariant quantity {H , H) = Hf/a diverges for 
r —f rn- Hence, the physical subspaces at the horizon 
and at infinity are one-dimensional, implying that black 
hole excitations do not exist in the generic case. 

So far, we were not able to exclude non-generic solu- 
tions by rigorous means: The first problem is that the 
potential in Eq. (^5|) is not manifestly positive (although 
numerical investigation suggest that this is the case). 
Furthermore, the boundary term arising in the integral 
argument given in Appendix ^ does not vanish at points 
where w' = 0. It is, however, clear that the potential is 
positive if £ is big enough. In this case the integral argu- 
ment applies, at least for excitations of the background 
solutions with one node. 



3. Conclusion 

Since Hi parametrizes the non-circular part of the 
metric, we have shown that there exists no non-circular 
deformations in the odd-parity sector. This completes 
the classification of the stationary odd-parity excitations 
of the BK solitons and the corresponding non-Abelian 
black holes. The only physically admissible non-Abelian 
stationary odd-parity excitations of these configurations 
are the rotating, electrically charged solitons and the two- 
parameter family of black holes found in [Q . All modes 
lie in the electric part of the distinguished sector £ — 1. 



VIII. DYNAMICAL PERTURBATIONS 

Stationary perturbations need to be analyzed in or- 
der to find equilibrium solutions which are infinitesimally 
neighbored to known static configurations, or to establish 
local uniqueness results. The linear stability properties 
of static background solutions are, however, described 
by non- stationary perturbations. In order to study their 
dynamical behavior by means of spectral theory, it is 
necessary to cast the perturbation equations into a sys- 
tem of pulsation equations, that is, into a wave equation 
whose spatial part is (formally) self-adjoint. Using the 
static EYM soliton or black hole background, our task 
is, therefore, to write the perturbation equations ( pSJ ) - 
(pll) in the form 



9' ; 



u = 0. 



(96) 



where ^ is a self-adjoint operator, containing spatial 
derivatives up to second order. For perturbations of the 
Schwarzschild and RN black holes this was achieved in 
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Eqs. (|69|), ([71|), (|80|) and (|8l|). For perturbations of non- 
Abelian background configurations, however, one needs 
to proceed differently: 

For I — (i.e., for radial perturbations), the above task 
was achieved in |15|, where it was shown that the static, 
spherically symmetric BK solitons and EYM black holes 
have exactly n unstable radial modes in the odd-parity 
sector, n being the number of nodes of w. 

For ^ = 1, we will show below that the metric pertur- 
bations decouple, and that the perturbation equations 
can be cast into a wave equation for the remaining YM 
perturbations, where the operator A is symmetric and 
positive. This will establish the absence of unstable odd- 
parity modes in the sector i — 1. 

For (. > 1, we were not able to derive symmetric equa- 
tions in terms of the gauge-invariant amplitudes H , A, B 
and C . However, a system of hyperbolic equations can 
be obtained as follows: By virtue of Eqs. (pO|), ( p4[ ) and 
( [45| ) one can express the time derivatives of the electric 
components Hq , Aq and Bq in terms of the magnetic com- 
ponents Hi, Ai, Bi and C and their first spacial deriva- 
tives. Equations (jSTI), @, (El) and (Eel) then yield a 



hyperbolic system of the form 



d^ 



92 



dt'^ dp^ 



K- 



d_ 
dp 



where u comprises the magnetic components, u — 
{Hi ,Ai,Bi, C), and where the radial coordinate p is de- 
fined as in (|8^). Unfortunately, neither the first order 
derivatives nor the potential V are formally self-adjoint. 

In ||lj] we have argued that the gauge-invariant am- 
plitudes used in the present paper are not suited to de- 
scribe dynamical perturbations, an exception being vac- 
uum gravity or self-gravitating Abelian fields. In order 
to obtain a symmetric wave equation one needs to in- 
troduce amplitudes which are adapted to the staticity 
rather than the spherically symmetry of the background. 
In terms of these new, curvature-based amplitudes, the 
odd-parity pulsation equations can be cast into the de- 
sired form (9q), as we have shown in [H. 

In the remainder of this section we present the dis- 
tinguished cases £ — and i = 1. In these situations 
the perturbation equations can be written in the desired 
form, since gravity can be decoupled for i = 1, while 
for ^ = only YM perturbations are present in the odd- 
parity sector. 



A. The pulsation equation for I — 

For spherically symmetric perturbations the pulsation 
equation is obtained from Eq. (p7 



d^ (R^da) + 2w^a = 0, 

for the gaug e- invariant YM amplitude a = a — A[v/w); 
seeEq. (E8). This equation is not regular at points where 



w vanishes. (The case where w vanishes identically was 
discussed in Sect. VC. Introducing the regular one- form 
wdv, and defining the potential $ 



vo^a — vj^a + 



vdw - 
by the equation uP'a — Sd(w$), we find 



-A$ 



, dw dw , 



1 



[w^ 
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$ = 0, 



where we have also used the background YM equation 
(pl|). For a static background we may assume a time 
dependence of the form exp(iwi), which yields 



92 w''^ a 



$ = w $, 



(97) 



where a prime denotes the differentiation with respect to 
p, and now $ = $(p). In order to overcome the difficulty 
that the potential is singular at points where w vanishes, 
one may perform the following super-symmetric transfor- 
mation: First, the operator on the LHS can be factorized 
and written as Q^Q, with Q and Q^ according to 



w op 



u, 



Q' 



-w 



d_l_ 
dp w 



where u is subject to the differential equation 



One may then write Eq. (|97|), 
^P = (3$, which yields QQ^ - 



2aw^ 



Q^Q^ = 0^2$, in terms of 
= w^vl'. Since lu'^^ = Q^^, 



there is a one-to-one correspondence between $ and ^, 
provided that u ^ 0. Furthermore, \E' is normalizable 
if $ is normalizable, and vice-versa, since (^ , '5) — 
(g$ , Q$) = {Q^Q^ , $) = cj2($ , $). 
The equivalent problem, QQ^'i' — tu^'^, reads 



dp^ 



+ —{-iw''~\) + 2u' 
R^ 



* = UJ^^, 



(99) 



where now the potential is regular, provided that u is a 
regular solution of Eq. (p8|). Since the function 

f 
^0 — wexp / u{p)dp 



satisfies Q^^o = 0, it is a solution to Eq. (99) for cj = 0. 
The key observation in |1^ is that there exists a solution 
to Eq. ( |9^ ) such that u/w"^ and u' are regular and ^0 
is normalizable. Since the factor w causes ^0 to have 
exactly n nodes {n being the number of nodes of w), 
this establishes the fact that the transformed pulsation 
equation (B9h admits exactly n unstable modes. 

It remains to show that each unstable mode of Eq. 
(P9| ) can be realized by a regular choice of the original 
amplitudes a and v. In order to see this, we first note 
that for a; 7^ the inverse transformation becomes 
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implying that the gauge-invariant combination w'^a is 
regular. Finally, one adopts the temporal gauge, at = 0, 
with respect to which Eq. (B7h yields 



d_ 
di 



ao 



2a 



«;$, 



implying that ap is regular. Using ^ — Q^ 
w~^{w^y + u^, as well as the i-component of 
*d(w$) in the temporal gauge, gives 



w'^a 



■—I' = * w^. 

ot w 

This establishes the existence of exactly n unstable 
modes of the original perturbation equations, since u/w 
can be chosen to be regular, implying that v is regular. 



B. The pulsation equation for ^ = 1 

We now show that for £ = 1 the gravitational perturba- 
tions can be expressed in terms of the YM perturbations, 
which yields a pulsation equation for the YM amplitudes. 
The gravitational amplitude h enters the perturbation 
equations ( pq ) only via thecoordinate-invariant combi- 
nation Fh = d{R~'^h). The crucial observation is that 
the second plus 2w times the third minus 2G times the 
first equation in {p&j yields the conservation law 



dt 



1 
2G 



R^Fh + R^ {Fa + 2wFt. 



Recalling the definitions Fa = da— F^ and Fi, — db—wFh, 
we find after integrating the above equation 



Fh= f (dd + 2w db 



Co ~, 
— *1 

i?2 



where cq is a constant, and / denotes the background 
quantity / = {R'^/2G+l + 2w^)^^. Using this expression 
for Ffi in Eqs. ( pq ) yields the symmetric, inhomogeneous 
equation 



d+ 



Gd 



-co*d 



/ 
2wf 



(100) 



for the gauge- and coordinate-invariant YM amplitudes 
a and b. The 2x2 matrices G and F are symmetric and 
given in terms of the background quantities by 



G = 



/i?2 / 4w;2 + R^/G 



F = 2 



-f w 
-2w 



-Aw 

-2w 
w^ 



-Aw 
2i?VG 



1 



where F is positive definite for w'^ ^ 1. (The case w = 1 



was already discussed in Sect. VI B 



The one-forms a and b may be expanded with respect 
to Schwarzschild coordinates t and p. 



Edt + Bdp, 



where E and B represent the gauge-invariant electric and 
magnetic YM fields, respectively. Using this in Eq. ( |100| ) 
gives, for a static background. 



d 



dp [a 



1 



G E' -B 



FE 



Co 



d_ 
dp 



f 
2wf 



d_ 
"dt 



1 



G(e' -B 



FB = 0, 



where the dot and the prime denote differentiations with 
respect to t and p, respectively. In particular, for sta- 
tionary perturbations, E = B = 0, we recover the facts 
that the electric and the magnetic perturbations decou- 
ple, and that B vanishes. 

For dynamical perturbations a homogeneous pulsation 
equation of the desired form is obtained as follows: Dif- 
ferentiating the first equation with respect to t and the 
second one with respect to p yields the relation 

FE = (FB)', 

where we have also taken advantage of the fact that the 
background is static. Using this to eliminate E from 
the second equation, we obtain the following two-channel 
wave equation with formally self-adjoint spatial part: 






w^ + 2 -3w 
-3w 2^2 + 1 



AG .^ 



wy 



1 




QB = 0, 



(101) 



where Q satisfies F = 2Q 

Q = 



1 —w 
-w 1 



Since the operator is symmetric and positive, we con- 
clude that the spherically symmetric EYM solitons and 
black holes have no unstable odd-parity excitations in the 
sector £ = 1. 



APPENDIX A: LINEARIZED RICCI AND 
EINSTEIN TENSORS 

In this Appendix we give the expressions for the lin- 
earized Christoffel symbols and the Ricci and Einstein 
tensors. As we have argued in Sect. II C| , it is sufficient to 
compute the perturbations in the ODG. The Christoffel 
symbols for an arbitrary (not necessarily static) spheri- 
cally symmetric spacetime are 



r 



be 



BC 



0, n, 



= 0, 

ra T^a 



be 



= 6^^R-^VcR, 
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where V denotes the covariant derivative operator with 
respect to the two-dimensional metric g defined in Eq. 
(|l|). In the ODG the metric perturbations ^ and their 
inverse become 

Sgab^5g''''^5gAB = 5g^'' = 0, 
SgAb = hSA, Sg^' = -h''S'', 

where all indices are raised with the background metric, 
i.e., h'' = g''''ha and S^ = g^^Se = R-^g^^SB- Us- 
ing this, and the background metric (p]), the perturbed 



Christoffel symbols, (5r'^^ = 
^9ai3;v), become in the ODG 



hg'^''{Sgai,-i3 + Sg0„;a 



Eventually, 5Gab = SRab — jR gabfj 5Rab in the 
ODG, which vanishes by virtue of Eqs. ( |Al| ) and the 
transversality of the spherical vector harmonics. 



5Gab = . 



(A5) 



The ODG expressions (|A3|)-(|A5|) together with Eqs. (g) 
and (^ evaluated in the ODG yield the desired formulae 

(ITil). 



APPENDIX B: STATIONARY SOLUTIONS OF 
RW-TYPE EQUATIONS 



5T% 



be 



SV^sc - h^^ 



0, 5T\c^S^9BcRh''VaR, 

^^ bc^ S V{bhc}, 



{BJC}i 



and 



"^ Be 



h,R-'g^^\/[BSD]. 



where V[bft.a] = \{^bha~^ahb) and V{f,ft.a} = \{Vbha + 
Vahh). It is now a straightforward task to compute 
the perturbed Ricci tensor in the ODG. Using SR^p = 






i5r'' n, one finds 



5Rab = , SRab = V'^haV^ASB} 



(Al) 



We discuss the conditions under which the stationary 
RW type differential equation. 



~drN[r)dr + ^U(r) 



* = 0, 



(Bl) 



admits only the trivial solution. The potential V{r) and 
the function N{r) are assumed to be non-negative for 
r/f < r < oo, and to have analytical expansions of the 
form 

N{r) - N^{r-rH) + 0{r-rH)\ V{r) = VH+Oir-m) 

in the vicinity of the horizon rn , and 

N{r) = 1 + O(r-i), V{r) = £{£ + 1) + 0{r-^) 



and 
SRAb^^^^ 



i?4V[b {ha]R- 



^(^RAR+{dR,dR))hb 



(A2) 



In order to simplify the expression for 6RAb we take ad- 
vantage of the background equations (0) to write RAR + 
{dR, dR) = 1 — ^g^^ Rab- Also using the transversal- 
ity of the spherical vector harmonics, g^^"^ aSb = 0, we 
find 



2g^^VcV[^5B] 



AS] =i{t+l)SA. 

A 



which we apply in the second term of Eq. ([A^). Fi- 
nally using the fact that 5GAb = 5RAb — ^SAhbg'^^Rfj.i, 
in the ODG, we obtain (with the identity R^g'^'^Rab = 
g^^GAB = R'^G%) the result 



SGAb = § ^ V"^ 



R^^[b{ha]R- 



-y(A + i?^G-) 



(A3) 



which holds in the ODG. The expression for SGab follows 
from the fact that in the ODG 6Gab — SRab, 



SGab = V^^^sjV^/ib 



(A4) 



as r — > oo, where i > 0, Ni ^ 0, V ^ Q. [In particular, 
the RW equation ( [l7|) meets the above conditions, and 
so does the Zerilli equation, describing vacuum perturba- 
tions with even parit y.] Under the above conditions the 
differential equation (Bl) has regular singular points |2l| ] 
at r — rn and r = oo, implying that ^ behaves like 



* = 



and 



* 



Piir-rn) 
log(r-r/f)(5i(r- rn) 

r-^-iP2(r-i^ 



for 



rn, 



^ir-') 



for r 



oo. 



Here, the Pi, 2 (a;) and Qi,2{x) are locally convergent 
power series with -Pi, 2(0) ^ and Qi,2(0) ^ 0. 

For non-negative N(r) and V{r) the standard integral 
argument. 



0< 



'■1 

r2 






1 



y*2 ^j. 



V^] 1'dr+[7V^a^*];j, 



implies that Eq. (Bl) has only the trivial solution ^ = 0, 
provided that the boundary term vanishes in the limit 
Ti -^ fH and r2 — > 00. In particular, this is the case if 
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the asymptotic flatness and regularity conditions imply 
that the solutions with Qi and Q2 must be excluded. 

As an example, stationary solutions of the RW equa- 
tion ( [17|) with £ > 2 can be excluded as follows: The vari- 
ation of the curvature components 5R™]^^^Ax'^ A Ax'^ — 

2R^d{H/R^)S/\^BSA] must be bounded, implying that 
^/i? must remain bounded as well. Hence, the solution 
with Qi is not admissible, and neither is the one with 
Q2, unless £ = I. 



APPENDIX C: EINSTEIN- YANG-MILLS 
BACKGROUND SOLUTIONS 

In this Appendix we recall the behavior of the static, 
spherically symmetric soliton and black hole solutions to 
the EYM equations (p4[)-(p6[) at the singular points (see, 
e.g., p3]): In the vicinity of the origin one has (with 
G=l) 



N{r) ==l-4&V2+0(r'*), 
5(r) = 5*0 [1 + 4 6^2 +0(r4)] 
w{r) = 1-br^ + 0{r'^), 



(CI) 



with parameters b — — ^^"(O) and Sq > 0. In the asymp- 
totic regime one finds 



N{r) = 1 - 2Mr-^ + Oir-"^), 
5(r) = 1 + 0(r-4), 
w{r) = ± [1 ~ 72Mr-i + Oir-"^)] 



(C2) 



with parameters M and 7. Finally, in the vicinity of the 
horizon the behavior is given by 



H , 



N{r)^^{r~rH) + 0{r-rHy 



S{r) 
w(r) 



■'H 



G 



l + —^{r^rH) + 0{r-rHf 
Gh , 



(C3) 



WH + -p^ir - rn) + 0{r - th) 



F, 



H 



where Fh — 1 — (w|f — l)^/r|j and G// — WH{w'jj — l)/rH- 
Here, the free parameters are wh = w{rH) and Sh- 



APPENDIX D: SU(2)-VALUED HARMONIC 
ONE-FORMS 

We construct a basis of su(2)-valued spherical har- 
monic one-forms which transform canonically under the 
angular momentum operator J, defined by 

JxT = CxT, 

where T is a tensor field over the spherically symmetric 
(pseudo-)Riemannian manifold {M,g). Here ^x denotes 
the Lie derivative with respect to an infinitesimal rotation 



X on Af . [In particular, for infinitesimal rotations in IR 
about the x'^'-axis, we define Jk = Jx^i where {Xk)rs = 

Using the commutator relations 

[Jx,d]=0, [Jx,*]=0, [Jx,dr]=0, (Dl) 

where X is an infinitesimal rotation, and hence a Killing 
field for g, it is not difficult to see that 

ydr, dY, -idr 

form a basis of spherical harmonic one-forms with total 
angular momentum £, where Y = Y^™' are the standard 
scalar spherical harmonics. The dual basis, C_i = Ye^,, 
C_2 = g^^'^sYej^, G3 = f]^^\7 bYcji^ is a linear combi- 
nation of the standard vector harmonics (see, e.g. p4|). 
(Here and in the following, e^, denote the standard basis 
fields of IR^ , e^ is the radial unit vector, and e^ is a basis 
of S"^ with dual basis 6 . The antisymmetric tensor tjab 
is defined by W^ = -q^O^ .) Since the operator d is par- 
ity preserving, while the operator * is parity reversing, 
C_i and C_2 have even parity, while Gg = g SaCb ^^^^ 
odd parity. Here Sa = "Hab^^Y denote the transverse 
spherical vector harmonics, g^^V bSa = 0. 

In order to construct su(2)-valued spherical harmon- 
ics, we use the isometry e^ ^^ Tk to identify IR with 
su(2), where the standard inner product on M^ corre- 
sponds to the normalized inner product Tr = — 2trace 
on su(2). Vector- valued tensors are identified with su(2)- 
valued tensors, and the operator d is defined by the ex- 
terior derivative V for vector- valued forms a = a'^v_^. 
da — v_{Da^ — v_^{doi'' 4- w* A a^), where cj* is the Rie- 
mannian connection with respect to the standard metric 
on H^ . (With respect to the standard basis, y_^ = e,, 
one has w*- — 0, and thus Da = da, whereas, with re- 
spect to the basis vectors e^ and e^, one finds lo^ — 6 , 
uj^j^ = —gABO^, and uj"^ — w^O The basis of su(2)- 
valued spherical harmonics becomes 

X^=YTr, X2=g^''TAVBY, X^=fl^''TAVBY, 

(D2) 

that is, Xj = G • r. Since the parity operator does not 
act on the inner index, Xi and X2 have odd parity, while 
X3 has even parity. 

A basis of su(2)-valued spherical harmonic one-forms 
is now obtained by the same procedure as above: Us- 
ing the commutator relations (Dl), with d general- 
ized as above, one obtains the nine basis vectors dX^, 
*dXk^ Xkdr. The decomposition da — TiDa^ — Va — 
TrgAB^^ A a^ of the total exterior derivative of a vector 
valued form a tangential to S'^ now yields the identities 



dXi ^ YdTr + TrdY , 

dAa ^ VX2 - TrdY , 
dXs ^VXs + TridY . 



(D3) 
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Furthermore, one has 



X = flXl + /2^2, 



iVXj, + VX2 + i{£ + l)YdTr = 0. 



(D4) 



By virtue of these identities one may also use the one- 
forms YdTr, TrdY, VX2 instead of dXi, dX2, idX^, 
or the one-forms YidTr, Tr*dY, *WX2 instead of *dXi, 
*dX2, dXj. In fact, the new sets turn out to be more 
convenient in order to derive the perturbation equations. 
In conclusion, the su(2)-valued spherical harmonic ba- 
sis one-forms with odd parity are 

Xidr, X2dr, Fdr^, r^dF, VX2, (D5) 

while the even parity basis one-forms are 

Xsdr, YidTr, TridY, *VX2. (D6) 

This is, however, only true for i > 1. For £ — 1 and 
£ = Q the above fields are not linearly independent. For 
£ = I the dimensions of both the odd and the even 
parity sectors are reduced by one, since V^Vsy^^"^' 
= -9abY'^'^=^^ implies VX2 = g^^VAVsYTcO^ = 
—YTAd"^ = —YdTr- For £ = 0, Y is constant, and hence 
X2, X^, and dY vanish. Specially, in the even parity case 
only idTr survives, which yields the spherically symmet- 
ric magnetic Witten ansatz for the gauge potential. 

It is also worthwhile noticing that the odd-parity ex- 
pansion (ra) of the metric perturbations can be obtained 
by "lo weri ng the inner index" and symmetrizing the one- 
forms ( |D6| ): 

Xsdr = g'^^SATBdr -^ Sg = SA^dr (g) 6"^ + 0"^ (g) dr), 
YidTr = YTAvJie^ -^6g = 0, 
TridY = TrSAO^ -^ Sg ^ SA{dr ® 9^ + 6^ ® dr), 
*VX2 = g^^TB^cSAO"^ -^Sg = ^{aSb]0"^ ® 0^ ■ 

In a similar manner the even-parity metric expansion can 
be obtained from the (odd-parity) one- forms (Dq). 



APPENDIX E: INVARIANT YANG-MILLS 
PERTURBATIONS 

In this Appendix we construct the gauge- and coordr 
nate-invariant amplitudes parameterizing the perturba- 
tions of the YM potential SA. Starting with Eqs. (27) 
^y and 



( P8| ) and (^9|), our aim is to show that the physical per- 
turbations for £ > Ij £ = 1 and £ — are given by the 
expressions (pO|), ( pi] ) and (|33|), respectively. 
Under YM gauge transformations one has 

SA^ SA + Bx, 

where D is the gauge covariant derivative with respect 
to the background connection (|l8), and x denotes the 
su(2)-valued scalar field parameterizing the gauge free- 
dom. For odd parity perturbations x is given in terms of 
two functions on M, 



where Xi and X2 are the odd-parity scalar isospin har- 
monics defined in Eq. (D2). 



Now using the identities (D3) and (D4) one finds 
BXi = TrdY + wYdTr, DX2 = VX2~wTrdY, the ampli- 
tudes defined in Eqs. (^) and ( p8| ) arc found to behave 
as follows under gauge transformations: 



a ^ a + d/i 
/? -- /3 + d/2 

M ^ M + /l - f2W 



for £> 1, 



and 



v + fiw-f2 for^ = l, 



V 

a 



V - 
a 



fiw 
h 



for £>l. 



(El) 



(E2) 



(E3) 



For £ — 1, one can introduce two gauge- invariant one- 
forms a and h, say, 



h = f3^ 









for 



1, 



(E4) 



which are well-defined unless the background configura- 
tion is the Schwarzschild black hole, w = I. The trans- 



formation laws (El) and (E2) imply that there exists a 
gauge for which the scalars fj, and v vanish. Moreover, the 
above definitions show that in this gauge the one-forms a 
and P coincide with the gauge invariant one-forms a and 
b. Since the perturbation equations are gauge-invariant, 
we may thus parametrize 5A^^^^' in terms of the two 
gauge-invariant one-forms a and b on M, 



SA^'- 



Xia + X2b. 



(E5) 



For £ > 1, we may proceed in a similar way and in- 
troduce two gauge-invariant one-forms and one gauge- 
invariant function as follows: 



a = a ~ d {^1 + wa) 
b = f3-da 
c=i>~'w{fi + wa) 



for £>l. 



(E6) 



It is again obvious from Eqs. (El) and (E3) that there 
exists a gauge for which /i and a vanish, and that the re- 
maining amplitudes a, /3 and v coincide with the gauge- 
invariant quantities a, b and c in this gauge. Hence, with- 
out loss of generality, we may set 



j^(^>i) = Xia + X26 + cYdr^, 



(E7) 



and consider a, b and c as gauge-invariant amplitudes. 

For .^ = 0, 5 A is parametrized in terms of the one- 
form a and the function u^ which transform according 
to a ^ a -I- d/i and v ^ v + fiw, respectively. The 
amplitudes combine into a gauge-invariant one-form 
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for e^o. 



(E8) 



where a coincides with a in the gauge for which ly van- 
ishes. (This gauge does not exist for the RN background, 
since v is gauge-invariant for w{r) = 0.) In terms of a 
one has 



SA('=o) 



(E9) 



So far we have parametrized 6A in terms of gauge- 
invariant ampHtudes, or, more precisely, in terms of am- 
phtudes which coincide with gauge-invariant ampHtudes 
in a certain gauge. However, these quantities are not yet 
invariant under infinitesimal coordinate transformations 
on the background. As the linearized Einstein and YM 
equations are invariant under these transformations, they 
will involve only coordinate-invariant combinations of the 
above amplitudes. In order to find these combinations, 
it remains to study the behavior of the gauge-invariant 
amplitudes a, b and c under the transformation 

6A^SA + CxA, 



where A is the background connection given in Eq. (18), 
and Cx denotes the Lie derivative with respect to the 
infinitesimal vector field X^" = -JR-'^S^jJ^'^^VbY, de- 
. In terms of the coordinate freedom /, 



fined in Eq. 
one finds 

CxA={\-w) 



i?2 



TrdY + VX2 



X.A 



J_ 

i?2 



(The most efficient way to establish this is to write Cx = 
dix + *xd, and to use zxdfJ — —R~^fdY and ix*dTr = 
R-^fX2.) 

The transformation properties of the one- forms a, (3 
and the functions /z, v and a defined in Eq. (|2^) are 
now immediately obtained. [For t = 1 one has to replace 
VX2 by -YATr and to use Eq. (||) instead of Eq. (|^.] 
For i > 1, the gauge-invariant quantities (^) transform 
as follows under coordinate transformations generated by 
X: 



■ a 

6- 
c - 



dw 



for ^ > 1, 



(ElO) 



'(1- 



while the transformation laws for the quantities ( E4 ) be- 
come 



■a-d(^) 

6-U.d(;^) 



for 1=1. 



(Ell) 



(There exist no allowed coordinate transformations in the 
odd-parity sector if £ = 0.) For £ > \, one may eventually 
use the transformation property (ra) of the metric pertur- 
bation K, K ^ K -|- /, to introduce the following gauge 
and coordinate-invariant amplitudes: 



A 



i?2 



(1 



B 



i?2 



dw 



C = C+—rw{l-W^ 



H.k-R^i{^) 



(E12) 



where we have also recalled the definition (10) of the 



coordinate-invariant metric perturbation one-form H . In 
the ODG {n = 0) these gauge-and coordinate-invariant 
amplitudes coincide with the gauge-invariant amplitudes 
a, h, c, and h, which reduce to the original amplitudes a, 
/?, 7, and h in the ODSG (k = ^ = cr = 0). 

For £ = 1 the gauge- and coordinate-invariant YM am- 
plitu des a re obtained by comparing the transformation 
laws ( Ell ) with the transformation property (||) of the 
metric perturbation /i, h ^ h + R^d{R^^f). This yields 
the invariant quantities a and b, defined by 



h 
R?' 



R^- 



(E13) 



APPENDIX F: LINEARIZED FLUX INTEGRALS 

The Komar expressions for the local electric and mag- 
netic charges, the local mass and the local angular mo- 
mentum of a stationary spacetime are given by the fol- 
lowing fiux integrals over a sphere with radius R: 



Qe[R) - 

M{R) = 
J{R) 



1 

1 

167rG Js^ 



^F, Qm{R) 

* (dgt^ A dx^) 



1 

47r 



F 



Sr 



Using the expressions (30) and (pll) for the gravitational 
and the YM perturbations, the linearized flux integrals 
are found to be 



SQm{R) = SM{R) = 0, 

6Qe{R)^Sne^-TR^{~*FA 

SJiR)^6iiS^oR^iFH, 



r*FB) . 



(Fl) 



where eg 



(0,0,1), e , J (Tl,i,0), and F^, Fa and 
Fb are defined in Eq. (|52|). Here we have also used the 
orthogonality of the spherical harmonics Y^™ and the 
expansions e^ - Y^"'e„, and S^"^ - g{dY^"',dY^°). 
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